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ABSTRACT 

We explore the prospects for using future supernova observations to probe the dark 
energy. We focus on quintessence, an evolving scalar field that has been suggested 
as a candidate for the dark energy. After simulating the observations that would 
be expected from the proposed SuperNova / Acceleration Probe satellite (SNAP), we 
investigate two methods for extracting information about quintessence from such data. 
First, by expanding the quintessence equation of state as wq{z) — wq{0) — Q;ln(l + z), 
to fit the data, it is possible to reconstruct the quintessence potential for a wide 
range of smoothly varying potentials. Second, it will be possible, to test the basic 
properties of the dark energy by constraining the parameters fig, wq and a. We 
show that it may be possible, for example, to distinguish between quintessence and 
the cosmological constant in this way. Further, when supernova data are combined 
with other planned cosmological observations, the precision of reconstructions and 
parameter constraints is significantly improved, allowing a wider range of dark energy 
models to be distinguished. 

Key words: supernovae: general - large-scale structure of Universe - cosmology: 
miscellaneous. 



1 INTRODUCTION 

There is no w strong evidence, from observations of ty pe la 

and 



supernovae (Perlmutter et al. 1999; 



Riess et al. 



CMB anisotropy (de Bernardis et al. 2002 



ically flat and domin 



19981) 



Netterfield et al 



2002; Pryke et al. 2002), that we live in a universe that is 



other forms of dark energy in so far as these can be param- 
eterised by a simple equation of state. 

It is our goal to explore what could be learned about 
the dark energy from future observations of type la super- 
novae (SNIa), such as may be possible with the pr oposed Su- 



ated b}r a nearly homogeneous perNova / Acceleration Probe (SNAP) satellite ( [Levi et al 



geomet 

rnmpmlpnt with npgativp pressnrp — the AarV pnprgy — which 

is causing the cosmic expansion to accelerate. The existence 
of dark energy has recently been confirmed, independently of 
the supernova data, by combining t he latest galaxy cluster- 
ing data with CMB measurements ( ^ilfstathiou et al., 2002 
Wang, Tegmark & Zaldarriaga 2002 ) . The most obvious dark 



energy candidate is vacuum energy, represented by the cos- 
mological constant A, which has pressure pA = — Pa- Other 
dark energy candidates, include a network of topological 
defects and the much-discussed po ssibility of quintessence 
(Caldwell, Dave & Steinhardt 1995), a spatially inhomoge- 
neous, evolving component which is usually represented by 
a scalar field evolving in a poten tial (an idea first intro- 
duced hv [Pechlos Ratra (19Ss|)V Tn thi 



200C). This satellite aims to observe roughly 2000 super- 
novae a year for three years, with very precise magnitude 
measurements and negligible systematic errors, out to a red- 
shift of z = 1.7. Observations of this sort will permit a 
very precise measurement of the magnitude-redshift relation 
m{z) and hence of the distance-redshift relation r[z), which 
will probe the expansion history of the universe. 

Current supernova observations already put weak con- 
straints on the dark energy density Q.q and equation of 
state Wq. If the universe is assumed to be flat (i.e., Q,k = 
1 — JIa/ — = 0) and wq is assume d to be constant in time , 
we have Q,q >0.5 and wq< - 0.4 ( Perlmutter et al. 1999 ). 
These constraints improve considerably when they are com 



lis paper wo wil 



bined with in dependent constraints from the CMB (Efs- 



focus o 1 nuiritessenco. tVioiiP-h sorno of the results a.nnlv to 



tathiou 1999|) or large-scale structure (Perlmutter, Turner 



& White 199t ) — most importantly, such combined data sets 
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give WQ< - 0.6. 

Several studies have examined the improved parameter 
constraints that may be possible with SNAP. For example. 
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Weller fc Albrecht (200l] , ^002| ) have found that SNAP can 
constrain a constant wq to better than 10% accuracy, al- 
lowing some models to be distinguished from a cosmological 
constant {wa ~ —1)- In addition, many authors have exam- 



Quintessence is usually parameterized by its equation of 
state Wq = Pq/pq- Looking at equations (jl|) and (Q, we 
see immediately that wq > —1. At late times, the coordi- 
nate distance to redshift z is given by 



lution of WQ with redshift (|Astior 2001; Barger & Marfatia 


2001 


; ( 


loliath et al., 2001; putcrcr & Turner 2001 




Maor, 


Bruste 


1 & Steinhardt 2001 


; Weller & Albrecht 2001 


20021). 



r(z) 



dz' 



Ho [^Mil + z')^ + nkil + z'Y + nQeC(-')]i/2 



(3) 



There is a general consensus that SNAP data alone will not 
be able to distinguish an evolving equation of state. Our 
results are in broad agreement with these studies. 

It may also be possible to perform a direct reconstruc- 
tion of wq{z) and of t he quintessence potentia l V from 
supern ova obser v ations (iHuterer fc Turner 1999| ; |Chiba &j 
Nakam ira 2000). 3aini et al. (2000) have attempted to per- 
form reconstruction from the current supernova data and 
find, unsurprisingly, that nothing definitive can be learned 
at present. The theoretical studies that have been done have 
often encountered difficulty in accurately reconstructing the 
properties of a given model from simulated sup ernova data 



where C,{z) = [3 Jg^(l -I- WQ{z'))d\n{\ + z')]. In this paper 
we will assume a flat (Q.k ~ 0) cosmology. The coordinate 
distance (3) is of fundamental importance since it fixes the 
magnitude-redshift relation probed by distant SNIa. 

The evolution of quintessence models is governed by the 
equation of motion 



Q + 3HQ + V' = 0, 



(4) 



(Huterer fc Turner 2001; Weller 



: Albrecht 2002| ). This casts 
doubt ripuu thy leliabihty of lecousLiuctioii. One of oui main 
goals ir i the present study is to develop a reliable method 
for producing accurate reconstructions from supernova data 
that is applicable to a wide class of quintessence models. 
Having done this, we will also want to explore whether such 
reconstt-uctions will be uriefrll. Thti ilifithod wti iisd fof feCon- 
structic 



where V' = dV/dQ. One difficulty with quintessence is 
that, in most cases, the initial conditions of this equa- 
tion must be fine-tuned for quintessence to dominate only 
just at the present epoch, as is required by observations. 
It is possi ble to address this problem by u s ing tracker 
potentials ( Zlatev, Wang fc Steinhardt 1999 ; Steinhardt 



Wang & Zlatev 1999), which admit attractor-like solu- 
tions, thus allowing quintessence to exhibit identical be- 
haviour at late times for a wide range of initial conditions. 
The energy scale of these models still requires fi ne-tuning 
to be consistent with observations, however (see Vilenkin 



Ju demunsLiaLeb the puwei ul usuig sapuinuva i 
vations to constrain cosmological parameters, and it high- 
lights the usefulness of combining these observations with 
prior knowledge from other cosmological measurements. We 
therefore also undertake an exploration of the parameter 
constraints that will be possible with SNAP and its combi- 
nation with other experiments, expanding on earlier analyses 
and discussing carefully what can and cannot be learned in 
this way. 

This paper is arranged as follows, in Section ^ we sum- 
marize the cosmological effects of quintessence pertaining 
to supernova observations. In Section ^ we simulate SNAP- 
like data for a quintessence-dominated universe and discuss 
the problems inherent in producing reconstructions from 
such data. We discuss a cosmologically parameterized fit- 
ting function for r{z) in Section H, and we use it for re- 
construction in Section In Section ^ we show how more 
general questions about the nature of the dark energy may 
be addressed by constraining parameters with SNIa data, 
and we discuss the long-term prospects for this sort of study 
in Section |^ We draw conclusions about theprospects for 
observational tests of dark energy in Section 



, ubsHi- 



(2001 ) f or further discussion of these issues). The fc-essence 



models (Armendariz, Mukhanov & Steinhardt 2000, 2001) 



also address the fine-tuning of initial conditions by adding 
to the quintessence Lagrangian kinetic terms that require 
quintessence to become dominant only after matter dom- 
ination. Non-minimally coupled scalar fields offer another 
possible explanation of why quintessence domination might 
be related t o the epoch of radiation and matter equality (see 
Bean (2001) and references therein). These models introduce 
extra complications with little extra motivation, however, so 
to keep things simple we will concern ourselves with tracker 
models exclusively here. We note that tracker models tend to 
evolve very slowly at late times, making it particularly dif- 
ficult to detect evolution and reconstruct potental shapes. 
So by focusing on tracker models here, we should be mak- 
ing a conservative assessment of the prospects for testing 
quintessence models with supernova observations. 

Inverse power-law potentials V oc provide the sim- 
plest tracker quintessence models, but such models are in- 
consistent with the current observational constraint wa ^ — 



0.6 dPerlmutter, Turner fc White 1999| ; [Efstathiou 1999[ ) for 
P > 2. To keep things both simple and realistic, then, we 
shall use a potential V oc as our standard test po- 

tential in this study. [Recently it has been suggested that 
current obse rvations may actually rule out all inverse-power- 



2 T RACKFT? QTTTIVTFSSFIMCF, CDSMnT.nnV law m odels ( Corasanti fc Copeland 2002| ; Bean fc Melchiorri 



We consider a scalar field (the qumtessence field) Q, evolv- 
ing in a potential V{Q). The field will be nearly spatially 
smooth, and for the purposes of this study, we ignore any 
small inhomogeneities. The density and pressure of quint- 
essence are given by 



PQ = ^Q' + V{Q) 



V{Q). 



(1) 
(2) 



2002 ) , but they nevertheless provide a useful theoretical test 
model.] We shall see next how a supernova measurement of 
r{z) can be used to probe the properties of the dark energy. 



RECONSTRUCTING THE QUINTESSENCE 
POTENTIAL 



Various authors (^tarobi nskv 1998 



[Huterer fc Turner 199£ 
Nakamura fc Chib^T999| ) have derived equations for recon- 
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structing the quintessence potential V{Q) and equation of 
state wq{z) from r{z). The equation of state of the quint- 
essence component is given by 



1 + wq{z) = 



(1+^) 



S^lMjl + zf + 2{cPf/dz^)/c^{df/dzf 
^ nuil + z)3~{cdr/dz)-^ 

and the quintessence potential is given parametrically by 



300 



200 



z 

CO 



LolQiz)] = 



(dr/dz) 



+ 



(l + z) d^f/dz^ 



3 [df/dzY 



(6) 



dQ 

dz 



= ± 



df/dz 

1 {1 + z)d'^f/dz'^ 
4-K (dr/dz)^ 




^QM{l + zf 
on 



1/2 



(7) 



Here we follow Huterer & Turner (1999) in using the di- 
mensionless quantities r — Hor, Q — Q/Mpi, and u){Q) = 
V{Q)/pc„t = ViQMpi)/{3Hi/8T:G). Thus, in principle, a 
measurement of r{z) will allow us to reconstruct the basic 
properties of quintessence. 

To test reconstruction, we choose a tracker potential 
V{Q), evolve equation (^) to obtain wq{z) and Q,q, and com- 
pute the Hubble-constant-free luminosity distance T>l (z) = 
(1 -|- z)f{z) from equation (^). The magnitude-redshift rela- 
tion is then given by 



Figure 1. The expected redshift d istribution of su pernovae ob- 
servations from the SNAP satellite (Levi et al. 2000), augmented 



by ground-based low-re dshift observations. Thi s distribution is 
similar to that given in Huterer & Turner (2001). 



m{z) = X +51og©L(2), 



(8) 



where M — M — 5 log Ho + 25 is the Hubble-constant-free 
absolute magnitud e, which we take to be = —3.45, the 
best-fit value from Efstathiou (1999). From this, we draw a 
Monte- 
points h 



It is important to note here a basic mathematical dif- 
ficulty with reconstruction. Looking at equation (Q), we see 
that reconstruction of V{Q) is impossible when the expres- 
sion in square brackets crosses zero. When this occurs at 
a given confidence level, we terminate the corresponding 
reconstruction contour for V{Q). Also, we note that the 
present-day value of the quintessence field Q is not fixed 
by equation (Q). But since this value has no cosmological 
effect, we may arbitrarily set Q = 1 at the present day. 
Various functional forms for f^^{z) have bee n recom 



'arlo sample of suporno-^ra magnitude redshift data 



SNAP JatolHtp ( [r.P^H Pt a1 9nnn[ ) 0,ir dmnlatoH Qampl 



200C 



mended in the literature for use in reconstruction (iHuterei 
fc Turner 2OO1I; IWeUer fc Albrecht 2002| ; |Chiba fc Nakamura 
2OOOD 



Saini et al. 



sists of roughly 2000 high-redshift supernovae from SNAP, 
in the range 0.1 < z < 1.7, plus 200 low-redshift points from 
ground-based searches, in the range < 2 < 0.2. The points 
have a mea surement error of 0.01 magnitude, plus an intrin- 
sic scatter (Efstathiou et al. 1999) of 0.018 magnitude after 
correction for the decline rate-luminosity correlation. The 
error in the redshift is negligible, and the redshi ft distribu- 
tion we use, ( Fig. |l|) , is similar to the one used by Huterer & 
Turner K^nni ) We neglect Sj^stematic errors in this analysis 



To show the difficulties inherent 
in reconstruction, we attempt reconstructions using two of 
them — a third order polynomial 



(2) 



3 2 

azz + a2Z ■ 



- a\z, 



and a Fade approximate 



2(1 -|- az) 



1 + &2 -f cz'^ ' 
both of wh ich have been suggested by 



(10) 



(11) 



Gi ven our set of N datapoints. we next choose a fitting 
function f^^{z). For a grid of points in the parameter space 
of r^^{z) we compute the likelihood function 



Huterer fc Turriei] 



(1999| , ^00 1| ). Reconstructions resuHing from fitting SNAP 



exp 



(mi 



2a? 



(9) 



where is the magnitude of datapoint i, 2i is its redshift, 
and (Ti is the error. The magnitude predicted by the fitting 
function is given by m^*(2) = M + 5 log O^^' (2). For the 
points in parameter space that fall within the 68%, 95% and 
99% confidence regions, we use f^'(2) to reconstruct wq{z) 
and V{Q) from equations (|^-0). We then draw confidence 
contours around the outermost reconstructed curves from 
each confidence region. 



like data with these fitting functions are shown in Fig ^ 
In applying the reconstruction equations here, we have as- 
sumed that SIm is known exactly. Evidently, the precision 
of the reconstructions will diminish if we allow for some un- 
certainty in Q.M- 

The reconstructions we obtain are inaccurate. For the 
polynomial fit, the actual value of wq lies outside the 99% 
confidence contour at low redshift, where the reconstruction 
is most precise, and the actual value of V{Q) lies outside 
the 68% contour at high values of Q. The Fade fit does 
significantly better, but the reconstruction of wq^z) is still 
slightly inaccurate at low 2. Reconstructions based on these 
fitting functions thus appear to be unreliable. Since our re- 
construction contours rule out the correct model, we cannot 
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Q/Mp, 



Figure 2. Reconstruction of V{Q) and wq{z) for a quintessence model with V oc Q~^. In the left-hand panels we h a.ve used the 
third-order polynomial (hol) as a fitting function, and in the right-hand panels we have used the Pade approximate (|ll[). The dashed 
lines — from innermost to outermost pairs — are 68%, 95% and 99% confidence contours, and the solid black lines are the actual model 
values. The heavy lines in the lower panels are lines of constant V, for reference. Inaccuracies in the derivatives of the fitting functions 
lead to inaccurate reconstructions of wq{z) and, for the polynomial fit, of V{Q). 



seriously claim that they rule out any specific model: they 
are too inaccurate to be useful. (Similar difficulties with the 
3rd-order polynomial fit were noted by Weller & Albrecht 



where a = dwo/dlna (Efstathiou 1999). Combining this 



(2002^n a preprint which appeared as this work was being 
completed.) 

This difficulty arises because the reconstruction equa- 
tions (^-Q) depend on the derivatives of r{z), rather than on 
r{z) itself, but the derivatives of the fitting functions do not 
necessarily resemble th e derivatives of r{z). (T his p roblem 
has als o been note d by Huterer & Turner (2001) and Weller 
& Albr jjcht (200^ .) The problem with the fitting functions 



suggested i n the literature is that, with the exception of 
the one in Saini et al. (200(j ), they are all arbitrary func- 



tions that have been chosen because they give good fits to 
r{z). However, we have no a priori reason to believe that 
their derivatives will resemble df/dz and d^f/dz^. Indeed, 
as Fig. ^ shows, the best fitting polynomial and Pade ap- 
proximate do not provide particularly accurate fits to the 
derivatives oir{z). Therefore we should not expect to obtain 
accurate reconstructions from these functions. To perform 
reconstructions reliably, we must first find a fitting function 
that we believe will provide not only a good fit to f{z), but 
also a good approximation to its derivatives. In essence, we 
want to fit our data with some physically motivated approx- 
imation to f{z). 



4 A PHYSICALLY MOTIVATED FITTING 
FUNCTION 

We can construct a physically motivated fitting function by 
noting that, for low z, the quintessence equation of state is 
often well approximated by 

w'q^ = WQ(ao) + aln(a/ao) = wq{z = 0) - aln(l -I- 2), (12) 



with equation (H), we arrive an at approximate expression 
for the coordinate distance, which we can use as a fitting 
function for reconstruction: 



dz 



[i + z'Y + ni^^{z')Y"' 



(13) 



where the approximate energy density contribution of the 
quintessence component is 

n;if*'(z) = nQ(l+2)='(^+™«(°» exp{-(3/2)Q[ln(l + z)f}.{lA) 

Also, by assumption, Q,m ~ 1 — (present-day values). 
Equation (^sj) is readily twice differentiable with respect to 
z; hence it can be used to reconstruct V{Q) and wq(z). A 
similar fitting function [with a linear, rather than logarith- 
mic, approximation to wq(z)] has been found to provide 
the best fits to r{z) of any o f the fitting functions pre viously 
considered in the literature (Weller & Albrecht 2001). Fig. H 
shows that equation ( |l^ ) also produces better fits to df/dz 
and d'^f/dz'^ than the polynomial or Pade approximate, sug- 
gesting that it may produce reliable reconstructions. This is 
entirely to be expected, since the fitting function is a di- 
rect approximation to the coordinate distance in terms of 
real physical parameters. Of course this function will pro- 
duce reliable reconstructions only to the extent that it is 
a good approximation to f{z), but we expect that it will 
approximate f{z) reasonably well, for a wide range of quint- 
essence mo dels, out to at l east the redshifts to be probed 
by SNAP ( [Efstathiou 1999[ ). It is worth noting that some 
models will not be well approximated by equation ( p^ ) — 
for example, t he pseudo-Nambu-Goldst one boson potentials 



considered by Weller & Albrecht (2002) (which can produce 



an oscillatory equation of state) would generally require a 
more complicated fitting function. Nevertheless, the analy- 
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Figure 3. Derivatives of the best fits to r{z) for various 
fitting functions, compared with the actual derivatives dr/dz 
and dP'rjdz^ (solid lines). The third-order polynomial fit [equa- 
tion ( [lO|. dotted lines], despite producing a good fit to r(g). gives 
an obvi 



jusly poor fit to the derivatives. It will therefore pro- 



duce poir reconstructions of V^fQ) and i/i(zV The Pade approx- 
imate fit [equation (|ll|), dot-dashed lines] provides more reason- 
able derivatives, but the small discrepancies between the fit and 
the actual derivatives still lead to inaccuracies in reconstruction. 
The fit based on a cosmological parameterization [equation jl^), 
dashed lines] provides the best fits to dr/dz and d'^r/dz^, as may 
be expected, since the fitting function here is an approximation 
to the actual coordinate distance. 



sis presented here is self-consistent: as long as equation (|13|) 
provides an acceptable fit to the data (reduced of order 
unity), we are unlikely to learn anything more by introduc- 
ing more complicated fitting functions. If there are signifi- 
cant residuals about the best fit, however, it may be fruitful 
to consider more complicated functions. 

It is also important to note that the use of equation ( p^ 
as a fitting function effectively reduces the problem of recon- 
struction to one of constraining the cosmological parameters 
Q.Q, wq{0) and a by likelihood analysis. Parameter estima- 
tion and reconstruction are separately interesting: the for- 
mer can tell us about the basic properties of dark energy in a 



model-iiudepeudeiiL way. 
specifickHy 



whereas ihn liiuer can tell uh mora 
Lial of a scdlai-field quiiiLessence. 



Since we are constraining some cosmological parameters in 
the reconstructions discussed here, it may be possible to 
combine supernova observations with other observations to 



better constrain the fit parameters and simultaneously to 
improve our reconstructions. For example, in a flat universe 
(which we assume), measurements of Hm from, e.g., galaxy 
clustering and the CIVIB fix fig via the constraint equa- 
tion Qq = 1 — Qm- If there are degeneracies between the 
parameters, then combining such measurements with super- 
nova data will serve to tighten the constraints on all the 
parameters, improving the precision of our reconstructions. 
Because we have used physical quantities to parameterize 
our fitting function for the coordinate distance, we can learn 
about r^^{z) indirectly from data other than the SNla obser- 
vations. In this way, a physically motivated fitting function 
allows reconstructions that are not only more accurate but 
also possibly more precise than reconstructions from more 
arbitrary fitting functions. 

Fig. ^ shows the sort of constraints on Q,q , wq (0) and a 
that might be expected from one year of SNAP data (2200 
SN). Clearly, strong and complicated degeneracies between 
the parameters make precise constraints impossible when we 
use the supernova data alone for estimation. It may be possi- 
ble to reduce these degeneracies significantly if we combine 
the supernova data with other measurements of the cos- 
mological parameters. For example, galaxy clustering mea- 
surements from the 2dF galaxy redshift survey, combined 
with measurements of CMB anisotropics, constrain to 
within roughly ±0.1 (2a err ors) in the case of a fiat uni- 
verse (Efstathiou ct al. 2002). To show the effect that such 
other measurements can have, we impose a Gaussian prior 
probability distribution on Q.m = 1 — ^q- Its effect on the 
contours, for various values of the standard deviation crn^.f , 
is shown in Fig. ^. Clearly the use of such a prior is signif- 
icant: if we impose (Jn^,j = 0.05, it is possible to constrain 
wq(Q) 7^—1 with 99% confidence, indicating that th e dark 
energy is not a cosmological constant. Other author s ( Asticr 
200l|; |Barger fc Marfatia 200l|; iGoliath et al. 200l|: iHutcrei 



Turner 2001 



Maor et al. 2001 



Weller & Albrecht 2001 



2002 ) have attempted to constrain a similar set of parame- 
ters using a SNAP-like dataset (using w\ = d-wg/dz where 
we have used a = dlnwQ fdlna). They also find that strong 
degeneracies between parameters make precise estimation 
impossible, unless prior knowledge of one of the parameters 
is assumed. Our results are in good agreement with these 
analyses.^ 

But is our Gaussian prior with crnj^^ = 0.05 realistic? 
The prior assumes that we will be able to measure flu to 
this accuracy independent of the equation of state, that is, 
that the measurement of Qm does not suffer from degenera- 
cies with Wq. In fact, measurements of flMh from the Sloan 
Digital Sky Survey (SDSS), combined with measurements 
of Q,Aih^ from the IMAP satellite's CMB measurements will 
be able to cons train Qm to nearly ±0.05, with no depen- 



dence on Wq (Eisenstein et al. 1998). IVIoreover, if flatness is 
assumed and this measurement is combined with measure- 
ments of the angular diameter distance to last scattering 
from the the IVIAP or Planck satellites, the uncertainty in 



Recently, ot her authors (Wang & Garnavich 2001; Wang & 



Lovelace 2001 ) have suggested, interestingly, that better infor- 
mation about dark energy may be obtained by constraining the 
function f{z) = e^(^' (where ((z) is as defined after equation (js])) 
at various values of z, rather than by attempting to parameterize 
and constrain w(z). 
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Figure 4. Likelihood contours for the cosmological parameters Hq, wq{0) and a, as defined in equation (tq). The contours axe calculated 
from a SNAP-like dataset for one year of observation, drawn from a quintessence model with V(Q) oc Q . The shaded regions indicate 
99% confidence regions, given various Gaussian prior probability distributions on Qm = 1 — f^Qi with standard deviation crfj^^. The 
contour lines inside each shaded region indicate the associated 95% and 68% confidence levels. Note that strong degeneracies between 
parameters make precise parameter estimation impossible in the absence of prior knowledge. In each graph, the two-dimensional likelihood 
contours shown are marginalized over the third parameter. 



Qm decreases to 0.01-0.03, with an additional strong con- 
s traint on wq a nd no degeneracy between the parameters 
( Hu et al. 1998 ). In addition, proposed X-ray or Sunyaev- 
Zeldovich effect surveys of galaxy clusters can constrain Qm 
to within an error o f 0.05 or better, with only a ve r y slight 



depend ence on wq (Haiman, Mohr fc Holder 2001 ; Weller 
Battye & Kneissl 2001). A proposed gravitational lensing 
survey could constrain flm to within 0.015 with minimal w q 



dependence (van Waerbecke, Bernardeau & Mellier 1999). 

The uncertainty an^.^ — 0.05 thus appears to be a con- 
servative estimate of the precision that may be obtained in 
the future. More optimistically, we could use a prior with 
c"njv/ = 0.015. Fig. ^ shows the constraints that are possi- 
ble in this case: with improved knowledge of Qm we can 
also obtain much tighter constraints on wq{0) and a. In no 
case is it possible to show that a 0, however, even if Qm 
is known exactly. Nevertheless, it is clear that combining 
SNAP parameter constraints with other measurements will 
be a powerful probe of the dark energy. We shall explore 
this in more detail below. 



5 RECONSTRUCTION FROM 

COSMOLOGICALLY MOTIVATED FITTING 
FUNCTIONS 

Prior knowledge of the parameters is also necessary if we are 
to use equation successfully for reconstruction. If we 
impose no prior probability distribution on our likelihood 



analysis (the "no prior" case in Fig. then reconstruc- 
tion of V{Q) fails for nearly all Q at all confidence levels, 
and wq{z) is only very weakly constrained. If, on the other 
hand, we impose a Gaussian prior with o"f2j^ = 0.05, we ob- 
tain the potentially useful reconstructions shown in Fig. ^ 
(Note that because we are explicitly constraining Qq — and 
hence Hm, which appears explicitly in the reconstruction 
equations — we have allowed Qm to vary over its allowed 
range in applying the reconstruction equations.) In partic- 
ular, V{Q) is constrained to be non-constant at the 68% 
confidence level. More importantly,the reconstructions are 
accurate: both the actual V{Q) and the actual wq{z) curves 
are contained within the 68% confidence contours. This is 
very promising: it should be possible, by combining SNAP 
data with expected measurements of Qm , to reconstruct the 
potential of an inverse-power-law quintessence model in an 
accurate and reliable way. 

Fig. ^ shows a reconstruction using the more optimistic 
prior = 0.015. In this case, we are able to constrain 

V(Q) to be non-constant at greater than 99% confidence 
with good accuracy, and we achieve impressive constraints 
on wq{z), although it is not possible to show that wq{z) is 
evolving. Because the reconstructions we obtain are accu- 
rate, it appears that equation is a better fitting func- 
tion for reconstruction than the arbitrary fitting functions 
we tried previously. But we cannot be entirely sure of this 
yet: when we attempted reconstructions with the polyno- 
mial and Pade approximate, we assumed exact knowledge 
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Figure 5. Reconstruction of V{Q) and wg{z) for a quintessence 
model with V oc Q~'^, using equation ( |l3| ) as a fitting function. 
A Gaussian prior is imposed on Sljvf, with standard deviation 
(JQm ~ 0-05 about the actual value. Both wq and V are accu- 
rately constrained, and we are able to constrain V{Q) const at 
68% confidence. 
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Figure 7. Similar to Fig. ^ except here we hold VLm fixed at 
its true value, rather than allowing it to vary. The fact that the 
reconstructions are accurate in this case indicates clearly that the 
cosmological fitting function ( [l3| ) is a better fitting function for 
reconstruction than the arbitrary ones we tried previously. 
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Figure 6. Similar to Fig. g, except here croj^ = 0.015. V(Q) is 
constrained to be non-constant at > 99% confidence and wq{^z) 
is constrained with good precision. 



of VIm. It is possible that our allowing to vary here 
might be obscuring inaccuracies. To test this, in Fig. ^ we 
perform a reconstruction from equation (^3|), in which we 
hold (and thus 0.q) constant at its true value through- 
out. The reconstructions are entirely acceptable: no obvious 
inaccuracies creep in when we fix fi^/. As hoped, the use of a 
physically motivated fitting function has given us a method 
for reconstructing the properties of quintessence in an accu- 
rate and reliable way. 



We note that 



Huterer fc Turner (2001 ) (HTOO) have pre 



respects from their analysis. First, we show here that equa- 
tion (^) produces good reconstructions of V(ff) as well as 
wq(^z). Second, whereas HTOO hold JIm fixed in their analy- 
sis, we show that good reconstructions are still possible when 
there is some uncertainty in Q,m- Third, HTOO use the Fisher 
matrix to draw their confidence contours, while we use an ex- 
act likelihood analysis. Most importantly, HTOO draw their 
simulated data from a model with wq exactly equal to Wq'^, 
which guarantees an accurate reconstruction (since the fit- 
ting function is the same as the original model). We show 
here that accurate reconstructions can be done even when 
is only an approximatio n to the actual equation of 



state. Weller & Albrecht (2002) have recently had some sue- 



viously attempted to constrain wq{z) using a fitting func- 
tion based on w'q^^z). Our work differs in several important 



cess in reconstructing wq{z) using a polynomial expansion 
of wq{z) to produce a fitting function, although they have 
also held fixed. Also, since the equation of state is not 
usually well approximated by a linear function of z, they 
must expand wq to second order to obtain accurate recon- 
structions, which requires one more fit parameter than our 
fitting function. We expect that this will lead to worse de- 
generacies in parameter space and that it may be difficult 
to obtain a useful reconstruction if Q,m is allowed to vary 
So it appears that, by fitting our data to equation (h 
we have found a way to reconstruct V(Q) and wq{z) in 
a reliable, accurate and fairly precise way. Recall, how- 
ever, that, in deriving equation (p^), our basic approxi- 
mation was to write the equation of state in the form 
w'q^{z) = wq{Q) — Q!ln(l -I- z). So what we are actually 
reconstructing when we use equation ( p^ ) as a fitting func- 
tion is the wq(z) that would he most likely to produce the 
data if wq{z) had the form of w'q^{z). Likewise, the V{Q) 
that we reconstruct is the potential that would have pro- 
duced this most likely w'q^{z). If we are careful to take these 
issues into account, it will be possible to make use of our re- 
constructions to distinguish between different quintessence 
models. As long as wq{z) is well approximated by w'^^{z) 
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over the redshifts we are measuring (i.e., as long as it is not 
dominated by nonlinear terms in ln[l + z]), we expect that 
our reconstructions will be accurate. Clearly models with 
rapidly varying wq{z) (e.g. , the pseudo-Nambu-Gold stone 
boson models discussed by Weller fc Albrecht (2002 )) will 
not be accurately reconstructed using equation (|l3|y7but in 
these cases we also would not expect an acceptable fit to the 
data. 

If we want to compare a given potential V{Q) to our 
reconstructions, we must first determine whether the wq{z) 
arising from it is well approximated by w'q'^{z). If it is not 
well approximated, we can say nothing useful about that 
particular potential using our reconstructions, since we can- 
not be sure that reconstructions of it would be accurate. If, 
on the other hand, the model in question is well approxi- 
mated by w'q'^{z), we can determine whether it is allowed 
by the data. First, we vary the parameters of the potential 
V{Q), including the present-day value of the field Qq, since 
this is not fixed by the reconstruction equations. If we do 
not find any parameterization that falls entirely within the 
reconstruction contours, then the model under considera- 
tion is disallowed by the data. If, on the other hand, we do 
find a satisfactory parameterization, we must then addition- 
ally check whether or not the resulting equation of state falls 
within the reconstruction contours for wq (z) for the redshift 
range of the data. If it does not, then the model is ruled out. 
Any model that passes both of these tests is consistent with 
the data. 

As an example, we test an inverse-power-law model with 
V oc Q~'^ against the contours shown in Fig. ^ (recall that 
these arose from a quintessence model with V oc Q~^)- We 
vary the parameters AI and Qo for the potential V{Q) = 
Af/IQ - (1 - Qo)f and identify a range of M and Qo for 
which the potential falls within the 68% confidence contours. 
This gives us a minimum and maximum allowed value for 
M. We then evolve the quintessence field in the potential 
V{Q) = M"^ 10^ , to get 'Wq{z) curves for both of these val- 
ues. (The transformation Q Q — (1 — Qo) has no cosmo- 
logical effect.) For inverse-power-law models wq decreases 
monotonically with increasing M; therefore, all allowed po- 
tentials will have wq(z) values that fall between these two 
limiting curves. As Fig. H shows, wq{z) falls outside the 99% 
confidence contours for all allowed values of M, and so mod- 
els with V oc Q~^ are ruled out by the data. The reason that 
the allowed potentials do not lead to allowed wq{z) curves 
here is that the calculation is not self-consistent. Each al- 
lowed parameterization of V{Q) includes a specific value for 
both M and Qo, but only M has any effect on the evolution 
of the quintessence field. For each allowed value of M, if the 
quintessence field does not evolve to the associated Qo at the 
present epoch, the resulting wqIz) may still be disallowed. 

Reconstructions using the fitting function equation ( p^ 
are thus clearly more useful than reconstructions using arbi- 
trary fitting functions. Using this equation, it is possible to 
determine, in a consistent way, whether or not a particular 
model is ruled out by the data — so long as w'q^{z) is a good 
approximation to wq[z) for that model. With the polyno- 
mial and Fade approximate fitting functions we tried before, 
we would have no reliable way of ruling out specific models, 
since we would have no good way of determining whether the 
fitting function accurately approximates the derivatives of f . 
Any apparent conclusion from arbitrary fitting functions like 
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Figure 8. Comparison of models with V{Q) = M"^ /Q^ to the 
confidence contorus in Fig. |6| In the lower panel, the dot-dashed 
line and solid line correspond, respectively, to the minimum and 
maximum values of M for which the potential falls between the 
68% confidence contours. The corresponding lines in the upper 
panel are wq{z) that arise from evolving quintessence in these 
potentials. Because wq decreases monotonically with increasing 
AI for inverse-power-law models, all allowed values of wq{z) fall 
between these two lines. Clearly, models with V oc are ruled 
out by the data. 



these could be spurious, the result of unaccounted-for inac- 
curacies in the fitting function or its derivatives. If we wish 
to produce useful reconstructions, it will be necessary to use 
a fitting function similar to equation (|l^). 

Testing individual models as described above could be 
useful, but because we lack any good theoretical motivation 
for model-building at present, testing the large numbers of 
possible models will be very time-consuming. For this rea- 
son, it may be more worthwhile to consider more general 
questions about the dark energy via likelihood analysis in 
parameter space, rather than to test models individually via 
reconstruction. 



6 PARAMETER ESTIMATION: 

DETERMINING THE PROPERTIES OF 
DARK ENERGY 

There are two especially important questions we would like 
to address concerning the properties of the dark energy. 
First, we would like to know whether the equation of state 
WQ differs from — 1 — that is, whether the dark energy is the 
cosmological constant or not — and second, we would like to 
know whether wq varies with redshift. In the previous sec- 
tion, we saw that it is possible to address these questions by 
performing a likelihood analysis over the parameters wq(0), 
a and Qq. In Fig. ^, we saw that it is possible, for an inverse- 
power- law quintessence model, to constrain wq{0) 7^ —1 at 
95% confidence using SNAP data alone. If we also include 
prior knowledge of fig , then it is possible to rule out A mod- 
els at even higher confidence. Even with exact knowledge of 
Qq, though, it is not possible to rule out a = for this 
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Figure 9. Likelihood contours for Qq and the effective con- 
stant equation of state from simulated SNAP data. The 
contours shown are for an inverse-powcr-law quintessence model 



with V{Q) ui. Q — ' ("PL model") and for a coHiiiolugical cousLaul, 
model (I 'A model"). Shaded contouro are 68%, 05% and 00% con 
fidence contours for one year of SNAP data (2200 SNe). The line 
contours define the same confidence regions for three years of data 
(6600 SNe). Clearly, it is possible to distinguish the quintessence 
model from a A model with high confidence for a full three-year 
SNAP dataset. 



model; it is not possible in this case to say whether wq is 
changing with redshif t. This result is in agreement with se v- 



eral previous studies (Astier 2001 

jHutorcr & Turner 2001 



Gohath et al. 2001 



Barger & Marfatia 2001 



Weller fc Albrecht 2001, 2002) 



Maor et al. 2001 



6.1 Two- parameter fits 

It is also clear from Fig. ^ that strong degeneracies between 
a and the other parameters lead to large uncertainties in 
our parameter estimation. Since we have not been able to 
say anything particularly useful about a, it might be more 
useful to leave a out of the analysis entirely for the sake of 
parameter estimation, focusing our attention on the value of 
Wq exclusively. To do this, we derive a fitting function for 
r{z) based on an effective constant equation of state w'q , 



-fit 
r 



dz' 



n3(H 



1/2- 



(15) 



Several other authors have considered this parameterization 
of the coordinate distance, and current supernova observa- 
tions constrain this parameter to be Wq^ < — 0.4 (Perlmut- 
ter et al. 1999). When th ese observations are comb ined with 
CMB (|Efstathiou 1999 ) or large-scale structure (Perlmut- 



ter et al. 1999) data 



the constraint becomes w'q < 



0.6. 



Recently there has been some suggestion that CMB data 



slightly favors w'^ > — 1 (Baccigalupi et al. 2002), although 



this conclusion assumes that some rather poorly constrained 
parameters are known exactly. 

Fig. ^ shows the likelihood contours that are obtained in 
the parameter space of equation ( |l5| ) for SNAP-like datasets 
drawn from a A model and a tracker model with potential 
V{Q) oc Q-^. Thefi gure shows the confidence contours that 
might be expected for one year and three years of SNAP 



data. No priors have been imposed on either of the param- 
eters. In the case of quintessence, with one year of data it 
is possible to constrain Wq^ —1 with 68% confidence, and 
with three years of data it is possible to do so with 99% 
confidence. Since wq = — 1 for a A model, a cosmological 
constant is ruled out at the same confidence levels. The fact 
that we can use SNAP to distinguish between a cosmolog- 
ical constant and quintessence is very encouraging. So far, 
however, we have only shown this for the specific case of an 
inverse-power-law quintessence model. It is not too surpris- 
ing that WQ can be distinguished from —1 in this case, since 
this model has present-day equation of state wq{0) ~ —0.6, 
relatively far from —1, when Qq « 0.7. It will be interesting 
to see if we can obtain similar results for a broader range of 
quintessence models — especially those with lower values of 

WQ. 

To test this, we consider two different tracker po tentials. 
The first is a potential inspired by supergravity (Brax & 



Martin 199£) 



ViQ) 



■ exp 



1 ( Q 

2 \Mpi 



(16) 



This potential exhibits tracker behaviour for P > 11. It 
also tends to produce lower wq than the inverse-power-law 
potential — for example for Qq ~ 0.7, we have wq (0) ~ —0.8. 
We also consider a quintessence model with a potential given 
by the sum of two power laws, 

where P2 = Pi x 10~®. Because the powers of the two 
terms differ by six orders of magnitude, this potential tends 
to produce values for the equation of state that are very 
close to —1 at the present epoch — when Qq ~ 0.7, we have 
wo(0) » —0.98. This t ype of potential was introduced by 



Steinhardt et al. (1999) to show that a very contrived po- 
tential is necessary for a tracker model to produce values of 
wq{0) very close to —1. We use this potential here to ex- 
plore the power of SNAP to distinguish quintessence from A 
models in the extreme case when wq{0) —> —1. Hereafter, we 
shall refer to simple inverse-power-law quintessence as a PL 
model, to the supergravity-inspired quintessence model as a 
SUGRA model, and to the sum-of-two-power-laws model as 
a 2PL model. 

Fig. ^ shows the likelihood contours that we obtain 
from SNAP-like datasets drawn from each of these models. 
In neither case is it possible to constrain Wq^ 7^—1; hence it 
is not possible to distinguish either of these models from a 
cosmological constant using SNAP data alone. (Our results 
for the A, PL, an d SUGRA models are in g ood a greement 
with the resuhs of [Weller fc Albrecht (2"00l| , |oo|).) Appar- 
ently SNAP data, on their own, can only be used to rule out 
A models in cases where wq differs significantly from —1. In 
section ^ though, we found that applying a prior probabil- 
ity distribution on Qq helped to reduce the errors on all the 
parameters of our fitting function. Using such prior informa- 
tion might help us to distinguish a larger set of quintessence 
models from a A model. We shall explore this possibility 
next. 

But first we note that there is a discrepancy between 
our quoted values of Qq and wo(0) and the likelihood 
contours shown in figures ^ and For example, for the 
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Figure 10. Similar to Fig. ^, except that here the UkcUhood 
contours are for data drawn from tracker models arising from the 
supergravity-inspired potential ["SUGRA" — equation ( |l6flj and 
the sum of two inverse power laws ["2PL" — equation (fTT])]. In 
neither case is it possible, using SNAP data alone, to distinguish 
the quintessence model in question from a A model. 



SUGRA model, we get present-day parameters Qq ~ 0.7 
and wq(0) « 0.8, but the likelihood analysis gives best- 
fitting parameters of fig ~ 0.6 and Wq^ ~ —0.95. Similarly, 
for the PL model, the actual parameters are Hq 0.7 and 
wq{0) ~ —0.6, but the best- fitting parameters are Qq ~ 0.6 
and = -0.7. 

It is important to keep in mind, though, that equa- 
tion ( p^ ) is nothing more than a fitting function to the co- 
ordinate distance. In this context, fig and Wq^ are merely 
fitting parameters, and our likelihood analysis simply finds 
the parameters that produce the best fit to the data. Be- 
cause our fitting function is only an approximation to r{z), 
the best-fitting parameters may turn out to be significantly 
different from the actual physical parameters they represent. 
(In the present example, the best fit is given by reducing 
both wq{0) and SIq from their actual values.) Nevertheless, 
these parameters could still be used to rule out a cosmolog- 
ical constant. In the case of a A model, the fitting function 
is an exact representation of f{z), so we would expect the 
best-fit parameters to match the actual parameters. Thus if 
we rule out wg = —1, we can say definitively that we are not 
dealing with a A model. But it is important to recognize that 
we can say nothing definitive about the actual values of the 
parameters, since the best-fit parameters do not necessarily 
reproduce them. 

Now we return to the idea of imposing prior probability 
distributions on our likelihood analysis. As in section |^, we 
assume that we have measured SI a/ = 1 — SIq with Gaussian 
error crn^j about its true value. In view of the issues just 
discussed, it is important to say clearly what we are con- 
straining here. The maximum likelihood parameters in this 
case are the parameters D,q and w'q that would be most 
Hkely to produce the SNAP data if wg were constant, given 
what we already know about Qg. Because we would still 
expect these most likely parameters to be the correct pa- 
rameters in the case of a A model, ruling out w'q = — 1 is 
still tantamount to ruling out the cosmological constant. As 
shown Figure ^ then, it is possible to rule out a A model 
with 99% confidence for the case of a PL quintessence model, 
by combining a year of SNAP data and a measurement of 




0.5 0.6 0.7 0.8 0.9 



Figure 11. Likelihood contours in Qg-Wg space for various 
quintessence models, assuming Gaussian prior probability distri- 
butions on f2]v/ = 1 — f^Q. The shaded contours indicate 68%, 
95% and 99% confidence levels assuming a prior with standard 
deviation crfiM ~ 0.05, and the line contours define the same 
confidence regions for crfijj = 0.015. All contours use a one-year 
SNAP dataset. Including such priors in our analysis allows us to 
distinguish the PL and SUGRA models from a cosmological con- 
stant, but they provide no discriminatory power in the case of the 
2PL model. 



Qm with (Tnjj = 0.05. In the case of the SUGRA model, we 
can rule out the cosmological constant at 99% confidence 
if we have ctqj^j = 0.015. This is encouraging; it appears 
that using SNAP in concert with other cosmological obser- 
vations will provide a strong test of the basic properties of 
dark energy. 

It is important to note, however, that the data cannot 
distinguish the 2PL quintessence model from a cosmolog- 
ical constant model, even with a very good measurement 
of Qm- Although this model is highly contrived, it is clear 
that SNAP will not be able to distinguish any arbitrary 
quintessence model from a A model using equation as 
a fitting function. As wg approaches —1, quintessence be- 
comes increasingly difficult to distinguish from a cosmolog- 
ical constant model, and our only hope for distinguishing 
quintessence from A is to probe the evolution of wg with 
redshift. 



6.2 Three-parameter fits 

To probe the evolution of wg with redshift, we return to 
the three-parameter fitting function that we used previously 
for reconstruction, equation ([l^), which was based on the 
approximate equation of state w'q^{z) = wg{Q) — aln(l -|- 
z). As we saw in Fig. ^, to get a tight constraint on a it 
was necessary to impose a Gaussian prior on Q.m with an 
optimistic value for (JUm of something like 0.015. Fig. 
shows the likelihood contours that would be expected for 
one and three years of SNAP data combined with such a 
prior. We consider the same three quintessence models as 
before. Only in the case of the SUGRA model is it possible to 
constrain a to be nonzero — and here only at 68% confidence, 
and only for a full, three-year dataset. 

It is hardly surprising that a 7^ is more readily shown 
for the SUGRA model. In this model, wg is evolving rela- 
tively rapidly at at late times, with a = —0.43, as opposed 
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Figure 12. Likelihood contours for tlie parameters Hq, wq and 
a, as defined in equation (jl3|. Tlie contours shown are for, from 
top to bottom, the PL, SUGRA and 2PL models. The shaded 
contours are for one year of SNAP data and the line contours 
are for three years of data, where in each case a Gaussian prior 
with cnjvf = 0.015 has been imposed on flj^.j = 1 — Oq. Only in 
the case of the SUGRA model is it possible to constrain a to be 
nonzero at any significant confidence level, and then only for a 
full, three-year dataset. 



to -0.11 for the PL model and -0.01 for the 2PL model. 
Clearly, SNAP will be much less powerful for distinguish- 
ing evolving from non-evolving dark energy models than it 
is for distinguishing quintessence from A models. In most 
cases, then, it appears that two-parameter fitting functions 
will be more useful than three-parameter fitting functions 
for constraining cosmology with supernova data. 



7 LONG-TERM PROSPECTS FOR PROBING 
QUINTESSENCE 

So far we have seen that SNAP will be a powerful tool for 
determining properties of the dark energy. But even with 
very good prior information, it may not be possible to dis- 
tinguish between two very similar-looking models, such as a 
2PL model and a A model. It has been shown that super- 
nova data do a better job of constraining some cosmological 
param eters when the su rvey is extended to higher redshifts, 
z > 2 ( lEfstathiou 1999| ). Ahhough SNAP wiU only track su- 
pernovae out to redshifts of 2 = 1.7, it could be possible in 
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Figure 13. Likelihood contours in Q.q-w'q space for the SUGRA 
quintessence model. Solid contours define 68%, 95% and 99% con- 
fidence levels for three years of SNAP data, and dashed lines de- 
fine the same regions when that data is supplemented by a sample 
of 30 supernovae at redshift 2 = 3. When the high-redshift sample 
is included, it is possible to rule out a A model {w'q = —1) with 
68% confidence. No prior knowledge of the parameters is assumed 
in either case. 



the long term to supplement SNAP data with higher-redshift 
data from powerful new telescopes like the Next Generation 
Space Telescope (NGST). To get some rough idea of the use- 
fulness of such observations, we simulate them by adding a 
sample of 30 supernovae at redshift 2 = 3 to our SNAP-like 
datasets. 

Fig. |l^ shows the improvement that can be made on 
a two-parameter likelihood analysis for the SUGRA model, 
when we supplement SNAP data with high-redshift obser- 
vations. Whereas previously it was impossible to rule out 
a A model using SNAP data alone, as shown in Fig. |l3[ 
when we supplement a three-year SNAP dataset with high- 
redshift supernovae, we can rule out a A model with 68% 
confidence, with no prior knowledge of the parameters. We 
see a similar shift in the confidence regions for the PL and 
2PL models, although the 2PL model still cannot be distin- 
guished from a A model. 

In the case of a three-parameter fit, the effect of adding 
a high-redshift sample is smaller, but it can still be impor- 
tant in some cases. Fig. |l^ shows the likelihood contours 
that would be expected in the case of the 2PL quintessence 
model, for a three-parameter fit to a three-year dataset, im- 
posing a Gaussian prior on flm with = 0.015. The high- 
redshift observations do not affect the confidence contours 
strongly enough to indicate a nonzero a, but they do con- 
strain wq{0) 7^ —1 at 68% confidence, ruling out a A model 
at that confidence level. Obtaining this result is somewhat 
surprising, since it was not possible to do so with a two- 
parameter fit. The reason it is possible here is that, since 
our three-parameter fitting function allows wq to vary with 
z, the parameters wq{0) and Qq are likely to be more accu- 
rately constrained, whereas the two-parameter fit required 
an offset in the best-fit parameters to make up for the z- 
dependence of wq. Three-parameter fitting functions, then, 
may be more useful than two-parameter fitting functions in 
some cases, even when they cannot be used to discern the 
evolution of wq. 

Finally, we note that Fig. ^ assumes a very optimistic 
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Figure 14. Likelihood contours for a three-parameter fit to 
SNAP-like data drawn from a 2PL quintessence model. Solid con- 
tours are for three years of SNAP data; dashed contours are for 
a three-year dataset supplemented with 30 supernovae at 2 = 3. 
With the high-redshift sample, it is possible to distinguish the 
2PL model from a A model. 



observational scenario. To distinguish the 2PL model from a 
A model, we must have three years of SNAP data with negli- 
gible systematic errors, plus thirty high-redshift supernovae 
from NGST, an independent measurement of Qm with ~ 5% 
errors, and extremely good evidence that the universe is flat. 
But although an observational programme that would lead 
to constraints like those in Fig. [l4| is ambitious, it is cer- 
tainly not impossible. The simulated SNAP dataset we used 
in this analysis has realistic statistical errors, and systemat 



ICS are e xpected to be negligible in comparison ( Levi et al^ 
200C()~}}GST, when launched, will almost certainly see at 
least some additional high-redshift supernovae. Also, as we 
have noted, our assumed prior constraints on ^a/ are re- 
alistic. Furthermore, the MAP and Planck satellites should 
provide very precise measurements of the cur vature of the 
universe [to within a percent or so for Planck ( [Tauber ct al 



200C |)]. [ The first steps have been taken toward observations 
that could lead to a measurement like Fig. Although this 
figure represents a best-case scenario, it also represents an 
achievable long-term goal. 



To answer our basic questions about the nature of the 
dark energy, however, cosmological parameter estimation 
may be more useful. If the dark energy equation of state 
is significantly different from —1, the SNAP data alone will 
be able to rule out the cosmological constant. Moreover, 
combining the SNAP constraints with independent measure- 
ments of Qm will allow us to rule out A models for all but 
the most extreme cases of quintessence, when wq is very 
close to —1. SNIa are less likely to provide evidence that 
WQ is evolving. Even combining SNAP data with very pre- 
cise measurements of 57 a/ will only provide evidence of an 
evolving wq in cases of extremely strong evolution. Three- 
parameter fits may still be useful, however: because they 
allow a more accurate determination of wq{0) than two- 
parameter fits, they may be able to rule out a cosmological 
constant in cases where a two-parameter fit cannot do so. 

The prospects are good for probing dark energy with 
SNIa data. SNAP, by itself, may be able to constrain cosmo- 
logical parameters with enough precision to rule out the cos- 
mological constant as a dark energy candidate. However, it is 
when SNIa data are combined with other observations that 
their true value becomes apparent. When combined with ex- 
pected future measurements of Qm, for example, SNAP data 
provide much stronger constraints on the quintessence equa- 
tion of state and may be able to produce useful reconstruc- 
tions of the equation of state and quintessence potential. In 
each of these cases, the combination of several observations 
tells us much more than any of the observations by itself, 
and in each case, the SNIa observations are an essential com- 
ponent. The measurements that may answer our questions 
about dark energy over the next decade or two constitute a 
difficult and impressive observational programme. Improved 
supernova observations are crucial to its success. 
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8 CONCLUSION 

In this study we have investigated the prospects for probing 
the dark energy, particularly quintessence, with observations 
of type la supernovae. We have examined tests for individ- 
ual quintessence models (reconstruction) as well as tests that 
address more general questions about the dark energy (pa- 
rameter estimation by likelihood analysis). In both cases, we 
find that data from the proposed SNAP satellite will pro- 
vide important information about dark energy, either alone 
or when combined with other cosmological observations. 

In the case of reconstruction it will be important to 
take extreme care in choosing a function to fit SNAP's mea- 
surement of the distance-redshift relation r{z). By using a 
fitting function based on a physical approximation to the 
coordinate distance, it is possible to produce accurate and 
reliable reconstructions from SNAP data combined with in- 
dependent cosmological observations. Such reconstructions 
might provide useful observational constraints as we con- 
struct models for the dark energy. 
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